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ABSTRACT 

The grand-design spiral galaxy M51 has long been a crucial target for theories of spiral structure. 
Studies of this iconic spiral can address the question of whether strong spiral structure is transient 
(e.g. interaction-driven) or long-lasting. As a clue to the origin of the structure in M51, we investigate 
evidence for radial variation in the spiral pattern speed using the radial Tremaine- Weinberg (TWR) 
method. We implement the method on CO observations tracing the ISM-dominant molecular com- 
ponent. Results from the method's numerical implementation-combined with regularization, which 
smooths intrinsically noisy solutions-indicate two distinct patterns speeds inside 4 kpc at our derived 
major axis PA=170°, both ending at corotation and both significantly higher than the conventionally 
adopted global value. Inspection of the rotation curve suggests that the pattern speed interior to 2 
kpc lacks an ILR, consistent with the leading structure seen in HST near-IR observations. We also 
find tentative evidence for a lower pattern speed between 4 and 5.3 kpc measured by extending the 
regularized zone. As with the original TW method, uncertainty in major axis position angle (PA) is 
the largest source of error in the calculation; in this study, where i5pA=±5°, a ~20% error is introduced 
to the parameters of the speeds at PA=170°. Accessory to this standard uncertainty, solutions with 
PA=175° (also admitted by the data) exhibit only one pattern speed inside 4 kpc, and we consider 
this circumstance under the semblance of a radially varying PA. 

Subject headings: galaxies: spiral - galaxies: kinematics and dynamics - galaxies: structure - methods: 
numerical 



1. INTRODUCTION 

The large angular size and clear spiral structure of the 
nearly face-on spiral M51 make it ideal for studies of the 
nature and origin of grand design spiral structure. Two 
scenarios dominate the discussion in the literature, each 
based on opposing theories: strong s piral structure a s 
a quasi-stationary density wave (e.g. iLin fc Shulll964f ). 
or as a transient feature du e to interac tion with nearby 
companion NGC 5195 (e.g. [riHwHIOTl . 

Observations of both the stellar and gaseous com- 
ponents reveal consistencies wit h the density wave 
interpretation at s ome level (se e lElmegreen et al.l 119891 : 
IVogel et all 119931: iRandl Il993f) In accord with the 
seminal study of iTullvl (|1974f ) which attributes the 
(transient) outer pattern to the interaction with its 
companion and proposes that the inner arms are likely 
spiral density waves a lso d riven by th e enco unter, 
lElmegreen et al.l (|1989f) and IVogel et al.l (|1993D find 
independent evidence for two different pattern speeds. 
If the strong spiral's corotation radius overlaps with the 



ILR of the o uter, material pattern, as is suspected (e.g 
iTullvl 11974 lElmegreen et al.lll989f) . this may indicate 
the stimulation of an inner spiral wavcmodc by the 
outer material spiral via mode-coupling (Ta gger et all 
Il987f ). But observationally it remains unclear whether, 
if driven by the outer pattern, the strong spiral structure 
is transient, or has surv ived a few rotation times (e.g. 

as speculated by IVogel "et al. 1991). 

In the simulations o f bothlSalo fc LaurikainenI (|2000ar ) 
and iHoward fc BvrdI (|1990f ) structure throughout the 
disk is well reproduced by multiple passages of the 
(bound) orbiting companion, and the nuclear structure, 
in particular, seems intimately related to the inward 
propagation of mul t iple ti dally-induced perturbations 
(ISalo fc LaurikainenI l2000bD . Pursuant to the study 
of iToomre fc Toomrd ()1972l) . such simulations of M51 
have proved indispensable for exploring and motivating 
scenarios in favor of short-lived waves. In the short-lived 
wave paradigm, propagating wave packets evolved from 
kinematic distortions in the outer disk may be swing- 
amplified (jToomrelll981l ). causing a strong response in 



the inn er disk. As predicted by ISalo &: Laurikainenl 
()2000b| ). the wave speed has a complex radial de- 
pendence, featuring a constant pattern speed for the 
dominant m=2 structure out to 1.2-1.8 kpc (depending 
on the disk mass assumed), followed by a superposition 
of structures described by a pattern speed that decreases 
with radius, down to '^10-20 km s^^kpc"^ by ^4.6 kpc. 
Knowledge of the pattern speed of the structure can, 
in principle, both distinguish between and reconcile 
the short- and long-lived wave scenarios, and so many 
studies have focuse d on measuring and characterizing 
this parameter (see l Elmeereen et a,l.lll989L lTullvlll974L 
Salo fc Laurikainenl l2000b' and iGarcia-Burillo et al.l 



19931 for example). The pattern speed of the outer 



spiral has long been proposed near 10-20 kms~^kpc~ . 
In the inner disk, application of the traditional, model- 
independent method of Tremaine & Weinberg (1984; 
hereafter TW) using CO observations yields a pattern 
speed r2p=38 kms~^kpc~ (jZimmer. Rand fc McGrawl 
|2004[ ). in general agreement with the determinations 
base d on r e sonan ce locations of lElmegreen et abl ()1989D 
and iTullvl (|1974f) (but higher than the pattern speed 
f]p=2 7 km s~^kpc~ found by iGarcia-Burillo et al.l 
Il993[ ). Although the TW analysis shows evidence for 
significant departures from the expected relation for a 
constant pattern speed in both the inner- and outer-most 
regions of the disk, the method cannot quantitatively 
account for any suspected radial variation of the pattern 
speed. 

The radial TW (hereaft er TWR) method 
(|Merrifield. R,and fc Meidtl [20061 : iMeidt et al.l [200ah 
should prove an invaluable resource in this regard, since 
with it we can characterize the angular speeds of distinct 
patterns and their possible radial variation. For the 
first time, we are able to observationally address issues 
related to the complex nature and persistence of spiral 
patterns and the connection, if any, between multiple 
pattern speeds in a single disk. 

Like its traditional counterpart, the TWR method, 
summarized in § 2, relies on the use of a kinematic 
tracer found to obey continuity. Here, we consider the 
ISM-dominant molecular component in the inner disk of 
M51 as traced by CO observations. In ijij l3.ll and l3.2l we 
describe these observations and review the arguments 
which establish their conformity with the assumptions of 
the method. In ij 13.31 we formulate the TWR quadrature 
and motivate the models developed for testing. 

Results of the regularized TWR calculation applied to 
the inner disk are presented in § 14.21 There we establish 
a best estimate for the pattern speed(s) of the bright 
spiral structure by considering the characteristics of 
solutions over a ±5° range of disk po sition angle s (PAs ) 
(§ ESI); according to the findings of iDebattistal (|2003[ ). 
we can expect uncertainty in the PA to be the dominant 
source of systematic error in the calculation. We also 
compare this estimate to other tested models of the 
radial dependence in § 14.41 

In an effort to authenticate the TWR estimate, in § 
14.51 we relate our measurements to other independent 
evidence for more than a single pattern speed in the 
inner disk of M51 and investigate the resonance locations 
and overlaps that they entail. We also consider our 
measurements in light of relevant findings throughout 



the liter ature, includi n g tho se of IShettv et al.l (|2007f ) (§ 
K6\i and iHenrv et all ()2003f ) (§ \AJ^. Final results are 
summarized in a conclusion section. 

2. THE TWR METHOD WITH REGULARIZATION 

The radial modification of the T\y m ethod 
(jMerrifield. Rand fc Meidtll2006HMeidt et al.ll2008L here- 
after MOB) delivers a derivation for radially-dependent 
pattern speeds measurable from observationally ac- 
cessible quantities. The so-called TWR calculation 
proceeds under assumptions parallel to those of the 
original method, namely that the disk of the galaxy 
is flat (unwarped); that the surface density of a disk 
component, which must obey continuity, becomes 
negligibly small at some radius and all azimuths within 
the map boundary (thereby critically yielding converged 
integrals; see below); and that the relation between the 
emission from this component and its surface density is 
linear, or if not, suspected deviations from linearity can 
be modeled. 

Departure from the traditional method (which as- 
sumes that the disk contains a single, well-defined 
rigidly rotating pattern) emerges by allowing that 
np=Up{r) - and the surface density of the tracer 
I](x,y,f) — 'S{r,(j) — rip(r)i). Integration of the conti- 
nuity equation obeyed by the tracer thereupon yields a 
Volterra integral equation of the first kind for flp{r), 
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(1) 



where x'{r,y) — ^r^ — y'^ (jMerrifield. Rand fc Meidtl 
|2006[) . This equation can be cast in terms of Xobs and 
yobsi the coordinates in the plane of the sky along the ma- 
jor and minor axes, respectively, and Vobs, the observed 
1.0. s. velocity, since for a galaxy projected onto the sky 
plane with inclination a^ x — Xobs, y — yobs/ cos a, and 
Vy = Vobs/ sin a. 

When the integral on the left of equation ([T]) is replaced 
with a discrete quadrature for different values oi y = yi 
and r = tj (represented in Figure 1 of MOB), equation 
III) takes the form of the matrix expression 

K.^n, = b, (2) 

with K an upper triangular N x N square matrix. This 
can be solved numerically for a total of two independent 
measures of flp{r), one from either side of the galaxy 
(y>0 and y<0). 

As described bv lMerrifield. Rand fc Meidtl (|2006D (and 
depicted in Figure 2 of MOB), solving Equation [2| by 
standard back-substitution results in the propagation 
of errors from large radii inward, whereby solutions in- 
escapably display noisy oscillations. As demonstrated 
there, applied first to Sb galaxy NGC 1068, this effect 
can be impeded most simply by adopting a relatively 
large bin width; the TWR solution in this case is found 
to decrease with radius, and yield a winding time esti- 
mate for the two-armed structure. 

But in general, as found in application to simulations 
(MOB), smaller radial bins are preferable to insure accu- 
rate assessment of radial variation. In addition, noisy 
behavior in solutions tends to be amplified when the 
quadrature extends out to the edge of the surface bright- 
ness (a requirement argued for by MOB), which not least 



imposes that the outermost bins generally cover the low- 
est S/N regions in the disk. Combined with a relatively 
small bin, numerical solutions as a result of inward er- 
ror propagation display a systematic offset in each bin 
between measurement and the actual value, preventing 
accurate determination of flp{r) (M08). 

As shown in M08, regularization provides an effective 
means of reducing the intrinsic propagation of noise in 
solutions while maintaining the precision required to ac- 
curately identify true radial variation. There, regularized 
TWR calculations were applied successfully to simulated 
disks featuring multiple pattern speeds in distinct radial 
zones as well as spiral winding. 

Following M08, then, we introduce a regularizing op- 
erator, or smoothing functional S, containing a priori 
informa tion in the manner of Tikhonov-Miller regular- 
ization (jTikhonov fc Arsenin|[l977l: lMilleiJll970l) into the 
X^ estimator minimized by solutions flj of equation ([2|) , 
whereby minimization returns smoothed solutions ac- 
cording to (in matrix form) 



(K^ ■ K + AS) • rj = K"^ • b 



(3) 



where the elements of K and b are Kij/ai and bi/ai, re- 
spectively (with errors ai representing the measurement 
error of the i*'' data point bi), and the parameter A con- 
trols the degree of smoothness achieved in solutions. De- 
tails for the full calculation and analysis can be found in 
M08; we proceed by highlighting only a few of the main 
precepts. 

By incorporating simple expectations from theory and 
observation into the smoothing S, the solution of cqua- 
tion[3]yields smoothed, testable models for flp{r). These 
models we restrict to simple forms and consider only 
polynomial solutions with constant, linear and quadratic 
radial dependence. (The elements of the smoothing S 
are associated with the minimization of the n*'' deriva- 
tive of ri(r) for each polynomial solution of order n.) 
These polynomial models can be incorporated into step- 
functions which parameterize the radial domains of mul- 
tiple pattern speeds (see M08). 

The best fit global solution constructed from the aver- 
age of like-model solutions from the two sides is estab- 
lished using the standard xt ix^ P^r degree of freedom) 
statistic, as in M08. (Note that an explicit assumption 
here is that all patterns in the disk are indeed global.) 
To summarize, once equation ([3]) is solved with a set 
of prescribed smoothings on each side, we use equation 
([2|) to generate a complete set of <v>i~ bi/{J T,dx)i for 
each global model. The xl fit of the model-reproduced 
to actual <v>i given global measurement error (t<"> 
(defined as the average of the individual errors (rf^^ for 
each slice) is then calculated for each. For this x^: we 
adopt the uniform weighting scheme advocated by M08. 

According to the M08 prescription where measurement 
errors for each slice refiect random noise in the data, for 
this analysis we assign errors (jf^^ that define the change 
in the measured <v> introduced by a change in the cho- 
sen flux cut-off in the first moment map. Specifically, the 
error in the <v> measured from a map with an na level 
cutoff is defined as 

^'^"^ = -^ [(< D >(„-l)a - <V >„<,)2 + (< 1. >(„+!)„ - <v >„„)' 

(4) 
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Fig. 1. — Zeroth (top) and first (bott om) moment maps of the 
M51 CO cube (originally presented in IShettv et al.l 120071 ). The 
horizontal bar in the top right corner indicates the physical scale. 
The y>0 quadrature generally covers the eastern half of the im- 
age (depending on the value of the PA adopted), while the y<0 
quadrature covers the western. 

for each slice i. The average of the individual errors then 
define the global measurement error across the entire disk 

where cr<"> = (E?fi c^^^''^)/2N (and N is the number of 

bins/slices used in the TWR calculation on a single side) . 
(Note that this error relates to the error at for each slice 
in equation ([3]) through {j T,dx)i.) 

Although these random errors are used in the 
goodness-of-fit criterion, the overall error in the measure- 
ment of Vtp{r) given by the best-fit global model solution 
is defined relative to systematic errors in the calcula- 
tion. Uncertainty in the assumed position angle (PA), 
for example, has the largest potential for introducing 
errors into <f>i, or conversely, the bi in equation ([3]), 
and is the d ominant source of error in TW calculations 
(|Debattista.2003; M08). We assess this error by testing 
the sensitivity of the solutions to departures from the 
nominal value for the PA (or inclination, for instance). 
Unless otherwise specified, in this paper all reported er- 
ror bars reflect the influence of PA uncertainty alone. As 
for inclination errors, apart from the change introduced 
in the pattern speed measurements through a change in 
sin a, these prove to be of little additional consequence 
to the accurate placement of radial bins defined in the 
quadrature (as suggested by M08) , despite the relatively 
low inclination (we adopt a=24°; see Table [Ij of the disk 
of M51. We therefore do not report this error, but in- 
stead note that a change in the inclination by i5q=±3° 
corresponds to a fractional variation of about 12% in the 
pattern speed estimates reported here. 

3. APPLICATION TO M51 

3.1. Observations 

In this paper, we consider the disk of M51 traced by 
high resolution C O observations. As described in its ini- 
tial pubhcation (jShettv et al.ll2007f ). the cube consists 
of the BIMA Survey of Nearby Galaxies (SONG) obser- 
vation together with several additional pointings which 
extend the map out to r~280" and provide higher angu- 
lar resolution in the central regions (see the beginning of 
§ 13.31) . A comp l ete de scription of the data can be found 
in She ttv et al.l ()2007[ ) . The 2-a zeroth and first moment 
mjips used in this analysis, derived from the full cube, 
I are shown in Figure [TJ Measurement errors given by un- 
certainty in the flux-cutoff are defined relative to maps 



at the 1- and 3-ct levels. 

3.2. Establishing Molecular Dominance 

The measured intensities and velocities in Figure [1] 
are suitable for use with the TWR method provided 
that the assumptions listed in the previous section are 
satisfied. While th e continuity requirement c an be 
parti cularly limiting , IZimrner. Rand &: McGrawl (|2004l ) 
and iRand fc Wallini (|2004( l argue that CO emission, 
the standard tracer of the molecular component of the 
ISM, suitably meets the TW assumptions for galaxies 
where the ISM is everywhere dominated by molecular 
gas. This is founded on the low true efhciency of star 
formation in spirals, which implies that only a small 
fraction of molecular gas is converted into stars on 
orbital timescales, while molecular dominance implies 
that the conversion of molecular hydrogen into other 
p hases of the ISM occurs at lo w leve ls . 

IZimmer. Rand fc McGra^ (|2004f ) applied the TW 
method under this premise, showing with CO and HI 
observations that the gas content of M51 is in fact 
dominated by molecular hydrogen where CO is detected. 
The CO observations used in this work can be similarly 
asserted to obey continuity: assuming a conversion 
factor between CO intensity and H2 column density 
X=2xl0^°cm~^[ K km s~^]~^, molecular hydrogen 
is found dominant over the majority of the CO cmit- 
ting disk (ro ughly i?<105"), where N(H9)/N(H1)-10 
(IShettvet al.l l2007). 

The possibility of variation in the CO-H2 
conversion factor has also been addressed by 
IZimmer. Rand fc McGrawl (j2004f ). who find in a series 
of tests applied to M51 that neither a linear relationship 
between metallicity and X-facto r nor arm-interarm 
variat ions at levels suggested by iGarcia-Burillo et al.l 
(|1993f ) produce a significant change in the derived 
pattern speed estimate. 

The negligible effect of radial dependence in X 
can be largely attributed to the cancellation of ax- 
isymmetry with TW integ r ation along each slice (see 
IZimmer. Rand fc McGrawl ()2004[ ) for a complete ac- 
count). Analogously, for the TWR method, as long as 
the metallicity changes negligibly over the width of a 
radial bin, we can expect little change in the results 
of calculations that assume an approximately constant 
X-factor throughout the disk. We have confirmed 
this to be th e case here; modeled according to the 
iBresolin et al.l (|2004[ ) metalhcity gradient 0.02 dex 
kpc~^, an increase in X with radius produces negligible 
change in the measured solutions. 

An arm-interarm contrast as suggested by 
iGarcia-Burillo et al.l (|T993), on the other hand, in 
general may not so readily translate from TW to TWR 
calculations inconsequentially. Currently, assessing 
the particular effect of azimuthal variation in the X 
factor on the results of the TWR method is beyond 
the scope of this work. Her e, we r ely on the results 
of IZimmer. Rand fc McGrawl (j2004f ) to assert that, 
to a first approximation, variation in X should not 
compromise the analysis as presented here. 

3.3. Defining the quadrature and developing testable 

models 



TABLE 1 
Parameters used in the TWR calculation. 



Parameter 



T^ 



Dynamical Center RA (a) (J2000) 

Dynamical Center DEC (6) (J2000) 

Distance 

Systemic Velocity (Vsys) 

Position Angle 

Inclination 



13'"29'"52^71 

47°11'42".80 

9.5 Mpc 

469 km s-i (LSR) 

170°±5° 

24°±3° 



Note. — The dynam ical center and inclination angle are 
adopted from lShettv et al] l|2007l) . Entries for Vsys and PA orig- 
inate from the tilted ring analysis of the first moment of the CO 
cube using the GIPSY task ROTCUR. 



In order to achieve as accurate a quadrature as possible 
and also limit errors caused by the misdesignation of any 
transitions (e.g. given the finite bin width to which so- 
lutions are confined; M08), we adopt a radial bin width 
Ar=0.23 kpc {D=9.5 Mpc). This corresponds to the 
limiting resolution (^4") of the map at the innermost 
radii. Since with the majority of our analysis of M51 we 
are most interested in characterizing the pattern speeds 
of the bright spiral structure, this choice is expected to 
yield high quality solutions for the pattern speeds in this 
region in particular. 

At the largest radii (and in interarm regions ) the reso- 
lution decreases to 6"-13" (|Shettv et al.ll2007D . Though 
in principle the quadrature can accommodate a non- 
uniform bin width, we maintain Ar=0.23 kpc through- 
out the disk and rely instead on the allocation of infor- 
mation administered by regularization. We assert that, 
even with our 4.5" radial bins, regularized TWR calcu- 
lations are prevented from oversampling the data as long 
as any distinct regions parameterized by the models are 
larger than the resolution. 

As assessed in M08, we can expect departures from the 
nominal values of the parameters appearing in Equation 
|2| to introduce non-negligible errors into the TWR solu- 
tions. Uncertainty in the major axis PA is the dominant 
source of systematic error in the calculation, resulting in 
errors on the order of 20% in TWR pattern speeds (M08). 
Since the kinematic parameters of M51 are notoriously 
difficult to constrain, perhaps the greatest challenge to 
the accuracy of our solutions lies in the the accuracy with 
which we can constrain the quadrature. 

To best equip the analysis in this capacity, then, we 
survey both our own derivations of the kinematic pa- 
rameters and those from the literature. For the coor- 
dinates of the center of rotation and the disk inclina- 
tion angle, for example, we rely on the values from the 
study of Shetty ct al. (2007). These wc then adopt in fits 
of a tilted ring model to the CO velocity field with the 
Groningen Image Processing System (GIPSY) program 
ROTCUR to determine the systemic velocity and the 
kinematic line of nodes (as well as the rotation velocity 
as a function of radius). The resulting parameters (listed 
in Tabic [1]) are consistent with most previous determi- 
nations. Note, however, that rather than adopting the 
range of PAs (170°to 180°) considered bv IShettv et~an 
(|2007D . we initially choose PA=170°±5°. This is princi- 
pally in order that our results are more easily compared 
with the majority of studies on M51, especially those 
which entail estimates for the pattern speed. Addition- 
ally, for this study we assume D=9.5 Mpc; with the al- 



ternative D=7.7 Mpc (more common to recent studies), 
all distances reported here decrease by a factor of '^0.2, 
while all pattern speeds reciprocally increase. 

According to the arguments in M08, we extend the 
unique quadrature established with the values in Ta- 
ble [1] out to the map boundary dzt/max in order to in- 
sure that all information critical for characterizing the 
patterns of interest is accounted for. Since the emis- 
sion extends (roughly East/West) out to ±ymax—^'^5'\ 
this defines the maximum radial extent of the quadra- 
ture Rmax = Vmax /cos a— 7.3 kpc, and hence the limit 
of integration along each slice ±Xi=^yRJr^ax — yi/ cos a. 
Though this does place Xi within the edge of the emis- 
sion at small \y\ (given the elongated emission in this 
map from N to S and the low disk inclination), the ra- 
dial range of the quadrature is still comfortably outside 
the radius where the integrals converge. 

As diagnosed by M08, with this fairly extensive 
quadrature solutions r2p(r) are at risk of regularization- 
induced bias. This bias is defined for the particular case 
when bins cover a region that displays only faint emis- 
sion, has little information from a strong pattern, or is 
suspected of sustaining multiple patterns; when these 
bins are prescribed an unrealistic model, the accuracy 
of the remainder of the solution can be jeopardized. 

Evidence for regions in the outer disk of M51 suscep- 
tible to regularization-induced bias are identifiable a pri- 
ori in the intensity map and its Fourier power spectrum. 
Later in § 14.21 where we address this bias and its signa- 
tures, we adopt the counter-measure developed by M08 
wherein the compromised bins are calculated without 
regularization. This imposes the additional parameteri- 
zation of a cut radius Vc on our model solutions, interior 
to which regularization proceeds as defined in cquation[3l 
(Sec M08 for a description of the calculation and analysis 
in this case.) In practice, we unregularize only as long 
as we can insure the sustained effectiveness of regulariza- 
tion in the rest of the calculation, given that an increased 
number of unregularized bins promotes the reintroduc- 
tion of unamendable propagating noise. 

In order to test for the possibility of multiple pattern 
speeds and/or winding, models additionally parameter- 
ize either single or multiple distinct radial domains over 
which the solution can vary as zeroth, first or second or- 
der polynomials. Though in general we test all possible 
models at each stage of the analysis, in some cases we 
restrict our consideration to only those polynomials for 
which the degree of freedom plus 3-4 bins does not exceed 
the number of bins in a given domain. 

4. M51: RESULTS 

4.1. Isolating the inner structure pattern speed 

As in previous applications of the method, we make use 
of a priori information to develop physically motivated 
models for ilp(r). And as with all such models, in order 
that they supply rigorous estimates we must also account 
for evidence suggesting susceptibility to regularization- 
induced bias. In the surface brightness (Figure [T]) and 
its Fourier decomposition (Figure [2) we identify a region 
outside r^A kpc, in particular, where both the surface 
brightness and power in the m—2 component are low, an 
indication that the information to be extracted there is 
potentially unreliable and difficult to constrain through 




Fig. 2. — Fourier power spectrum of the moment zero map shown 
in Figure [T] Modes up to m =4 are plotted as a function of radius 
with hnes for m =1 in black dash-dot, m =2 in red dash, m =3 in 
green solid, and m =4 in blue dash-dot-dot-dot. 

modeling. 

This is manifest in solutions for which regularization 
is employed throughout the full extent of the emission. 
Bins at large radii in the lowest-x^ solutions exhibit a sig- 
nificant degree of variation in their modeling, confirming 
that constraining the outer pattern speed is difficult. Ac- 
cording to the conclusion drawn by M08 in tests of the 
regularized TWR calculation on simulations, this chal- 
lenges the accuracy with which all inner bins can realize 
the true pattern speed. We therefore initially consider 
models which parameterize a cut radius rc=4.1 kpc, be- 
yond which all bins are calculated without regularization. 
In testing, we find this cut radius to coincide with a clear 
minimum in the x^, with all other best-fit parameters 
held fixed. 

A second, shallower minimum at rc=5.3 kpc is also 
compelling, and we consider its parameterization in mod- 
els of ilp{r) in the analysis that follows, as well. This 
location may well be reasonable for the separation of 
the patterns given that it seems to match expectations 
for the location where the outer, mat erial pattern be- 
gins. We cite in partic ular the study of 'Elmcg reen et ahl 
(|l989i) who, hkelTullyl (|T974 ). argue that OLR occurs at 
the termination of the bright, inner spiral structure and 
in the pretext of mode-coupling therein identify an over- 
lap between the CR of the inner pattern with the ILR of 
a 10-20 km s~^kpc~^ outer patter n (e.g. that f irst pro- 
posed for the material pattern bv iTullvl ()1974l )). This 
places the the innermost extent of the material pattern 
at r^ 6.0 kpc (adopted in to the distance convention used 
here). IVogel et al.l ()1993D also argue for a similar corota- 
tion radius based on observations of streaming motions 
in the ionized gas component of the ISM. And while this 
does not locate the inner extent of the outer pattern, 
it nevertheless implies that the outer arms are separate 
from, and have a lower pattern speed than, the inner 
arms (jVogel et al.l 11993"). C onsequently , it is consistent 
with the conclusion of Elmegreen et al.l (|1989f ). 

A transition from an inner to an outer pattern is 
also recognizable in the tida l perturbation-only model of 
ISalo fc: Laurikaineiil (|2000bl ). There, an independent spi- 
ral pattern with corotation near r=4.6 kpc is found to 
be followed by structure at the lower 10-20 km s^^kpc"^ 
pattern speed. Any resonance overlap, however, they ar- 
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Fig. 3. — The best-fit regularized solution with rc=4.1 kpc for 
PA=170°-|-/-5°. For this solution, bins exterior to r=4.1 kpc (not 
shown) have been calculated without regularization. Dashed red 
lines represent the difference from solutions derived with a two- 
pattern speed model at PA=165° and 175°. Horizontal error bars 
represent the dispersion in rt,i and rt,2 from PA to PA. The values 
in the zone of the bright spiral structure correspond to f2p i=90- 
27/+20 km s-^kpc"^ out to rt,i=2.1±0.3 kpc and np,2=50-|-9/-ll 
km s~^kpc~^ out to re- Curves for fl, Q±k/2 and n-K/4 (see § 
I4.5|l are shown in gray. 

gue is likely coincidental since the value of the higher, 
inner speed is associated with the maximum in f2-K/2, 
while the lower speed is determined mainly by external 
forcing. 

In addition, as revealed in the sections to follow, while 
results with rc=4.1 kpc indicate much higher speeds 
(-50-100 kms-^kpc"^) than the TW method {Qp=S8 
kms^^kpc^ ), with rc—5.i kpc solutions measure a 
much lower speed exterior to r~4 kpc, at least quali- 
tatively more consistent with the gross overall speed es- 
timated with the TW calculation. 

4.2. Best-fit models 

When we minimize the influence of the suspected ma- 
terial pattern in the TWR solutions by calculating the 
bins in the outer zone without regularization, we in prin- 
ciple maximize the leverage on the inner structure. In 
doing so, we find the data at the nominal PA to be well 
fit by two distinct pattern speeds interior to rc=4.1 kpc. 
The overall pattern speed solution with PA uncertainty 
SpA=±5°, to be discussed at length below, is represented 
in Figure [3l Following the treatment of M08 for con- 
structing errors on the measurement from a particular 
observational scenario, error bars represent the disper- 
sion of the parameters in the best-fit solution derived 
with a two-pattern speed model at PA=165°, 170° and 
175°. As will be discussed further in § 14. 5i these two 
pattern speeds both end at corotation, within the uncer- 
tainties. 

To quantify the relative benefit of the two-pattern 
speed solution, in Table [2] we list the xl estimate for 
several model solutions calculated over all slices in the 
TWR quadrature. In this table we also consider a xl 
over slices in the zones 0<|y|<2.3 kpc and 2.3<|y|<4.1 
kpc at the nominal PA. We expect the <v> for slices 
in each of these zones to predominantly reflect measure- 
ments in the radial bins r=\y\, so these separated xt 
(labeled Xu s hereafter) should provide a fair compari- 
son of il.p{r) from model to model at these radii. (Note, 



Fig. 4. — The best-fit regularized solution with rc=5.3 kpc 
for PA=170°-|-/-5°. For this solution, bins exterior to r=5.3 kpc 
(not shown) have been calculated without regularization. Dashed 
red lines represent the difference from solutions derived with a 
three-pattern speed model at PA=165° and 175°. Horizontal 
error bars represent the dispersion in rt^i and rt,2 from PA to 
PA. The values in the zone of the bright spiral structure cor- 
respond to r2p,i=96-26/+16 kms~^kpc~ out to rt^i=2.3it0.1 
kpc, r2p_2=51-|-7/-ll kms~^kpc~ out to rt, 2=3. 9it0. 4 kpc and 
rip,3=23-7/ +6 k ms~^kpc~ out to Tc- Curves for fl, Q,±k,/2 and 
Q,-k/A (see § 14.51 1 are shown in gray. 



however, that all outer bins also appear in the <v> re- 
produced by solutions in these zones.) 

The xl fit over all slices principally suggests that the 
two pattern speed solution and the single, constant speed 
solution yield better agreement with the data than the 
quadratic solution. For the former two solutions, the xt 
values are nearly indistinguishable at this PA (see the 
last column in Table [5]) . From the X^ g, on the other 
hand, it is clear that inside r~2 kpc the fit of the two 
pattern speed solution is significantly better than that of 
the constant speed. 

Comparisons with the <v> reproduced by the best-fit 
two-pattern speed model clearly demonstrate the incom- 
patibility of the constant speed model, as shown on the 
left side of Figure [5l There, the <v>i (top) and resid- 
uals (bottom) at each slice position reproduced by the 
best-fit solution are plotted along with those reproduced 
by the best polynomial solution with constant pattern 
speed 51p=55 kms~^kpc~ calculated over the same ra- 
dial zone. In the latter case, a greater departure from 
the measured values is readily apparent at slices inside 
Ij/I'^2 kpc as compared with the best-fit solution, which 
transitions from an outer speed of f2p=50 km s^^kpc^ 

to an inner speed of r2p=90 km s~^kpc~ at rt=2.3 kpc. 
In fact, the constant solution fits the data better than 
the two-speed solution in only 2 of the 23 such slices. 

In constrast, the quadratic solution with rc=4.1 kpc 
at PA=170°, which declines smoothly with radius from 
—95 kms^^kpc" , grants nearly comparable agreement 
with the measured values that the two pattern speeds 
entail, over a number of slices. However, the fit of this 
solution weakens at slices between |y|— 2.3-4 kpc (clear 
from the xt s)^ raising its x^ well above that of the two- 
speed model. 

By comparison, then, it would seem that the two- 
pattern speed solution presents the best fit for slices at 
both small and large radii (in slices |i/|<4 kpc). 
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Fig. 5. — (Top) Plots of solution-reproduced (dots and crosses) and actual (open circles) integrals <v>i=bi/ J Sdz as a function of slice 
position y at PA=170°. The values associated with the best-fit two-pattern speed solution calculated with rc=4.1 kpc (black dots) are 
plotted along with those of the order-zero polynomial solution with constant pattern speed Qp=55 kins~^kpc~^ (red crosses; left panel) 
and the second order polynomial solution (red crosses; right panel). (Bottom) Plots of the residuals in <v> reproduced by the solutions 
considered above (left: black dots for the best-fit, red crosses for the constant solution; right: black dots for the best-fit, red crosses for the 
quadratic solution). The adopted global error cr^"^ is shown in the upper right in each plot. Only those slices which show a contribution 
from bins inward of rc=4.1 kpc are shown. 
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3.64 
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Note. — Representation of the goodness of fit for model solutions with either rc=4.1 kpc (left) or rc=5.3 
kpc (right) at the nominal PA=170°. The x^ estimates in the zones 0<|j/|<2.3 kpc, 2.3<|j;|<4.0 kpc, and 
4.0<|j/|<5.3 kpc, labeled A, B, and C, respectively, are listed along with those for all slices, where all values 
are calculated as a reduced x^ difference of the model-reproduced to measured <v>i. The xt s i^i each 
radial zone is normalized by the number of bins in that zone minus the number of degrees of freedom, and 
so represent a goodness-of-fit distinct from a x^ over all slices with which we judge the whole solution. 
The multi-speed solutions with rc=4.1 and 5.3 kpc correspond to the best-fit two- and three-pattern speed 
solutions, respectively. Pattern speed estimates in units of kms~^kpc~ in each zone are also listed. For 
quadratic solutions, the values in the first and last radial bins are indicated. 



4.2.1. Extended models 

When the regularized zone 4.1<r<5.3 kpc is included 
in solutions, the best-fit solution once again measures 
two pattern speeds inside 4 kpc, but now a third, distinct 
pattern speed is also parameterized. Figure |4] plots this 
best-fit solution, where, again, error bars are defined by 
the dispersion in the lowest-x^ solutions derived with a 
three-pattern speed model at PA=170°, 165°, and 175°. 

According to the xt s i'^ Table [2] for solutions with 
rc=5.3 kpc, in both inner and outer zones this best-fit 
solution is superior to a two-speed solution which tran- 
sitions from a single constant 62 kms'^^kpc" pattern 
speed inside 3.2 kpc to a lower 20 kms^^kpc^ speed 
and also, once again, to either of the two polynomial 
solutions considered here. 

The quadratic solution (decreasing from 180 
kms^^kpc^ ) fits relatively well in the innermost 
bins, but overall the fit is now less comparable to that 
provided by solutions measuring three pattern speeds. 
The agreement between the constant model solution- 
reproduced <v> and the data also weakens, relative to 
the best-fit solution, especially inside r=4 kpc. This 
can be attributed to the decrease in the value measured 
with the constant model, from flp—55 kms^^kpc" 
with re— 4:.l kpc to ftp— 27 kms^^kpc^ with rc=5.3 
kpc. Incidentally, this is a clear indication that not only 
is the pattern speed in the zone 4<r<5 kpc lower than 
rtp—55 kms~^kpc~ , but as such undeniably influences 
all bin values calculated inward with this type of model, 
thereby interfering with accurate measurement interior. 

In contrast, when the zone 4<r<5.3 kpc is distinct 
and isolated, multi-speed models are nearly free of 
such inaccuracy. In both the two- and three-speed 
solutions with rc=5.3 kpc, the inner pattern speeds are 
nearly identical to the values measured in solutions with 
rc=4.1 kpc. This seems to suggest that, despite the 
evidence in Figures [1] and [2] that information beyond 
rci4 kpc is not conducive to modeling and extraction, 
the determination in the third zone is fairly accurate. 

This equivalence inside r~4 kpc to the pattern speeds 
measured in the solutions with rc=4.1 kpc derives in 
practice through the parameterization of a transition at 
r~4 kpc. This establishes an identical radial domain for 
the inner speeds in the solutions with rc=5.3 kpc and 
4.1 kpc. The transition rf_i^2.3 kpc in the three-speed 
solution as such yields the greatest similarity to the 
best-fit two-speed solution with rc=4.1 kpc. 

Moreover, according to Equation [21 rigorous measure- 
ment inside 4 kpc in principle also owes to accurate 
measurement for the pattern speed in the zone 4<r<5 
kpc. As inferred above, the measurement flp^3=20 
kms^^kpc^ is in fact lower than all measurements 
interior. Presumably, it is the value in this zone that 
contributes to the measur ement of the rather low 
TW value 38 kms^^kpc^^ (jZimmer. Rand fc McGrawl 
[200l . 

Even if the measurement for a distinct speed in the 
zone 4<r<5 kpc is a good description of the pattern 
there, since all x^ s-re lower in solutions with rc=4.1 kpc 
than with rc=5.3 kpc, we take this as an indication that 
calculating bins in the zone 4.1 <r<5.3 kpc without 
regularization does not reintroduce noise into solutions. 
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Fig. 6. — Best-fit regularized solutions at three PAs {PA=165°, 
170°, and 175°). Solutions calculated with both rc=4.1 kpc and 
5.3 kpc are plotted, with the latter shown in thinner line. In these 
solutions, bins exterior to re (not shown) are calculated without 
regularization. The solutions from PA=165°, 170°, and 175° are 
shown in blue dash-dot, green dash and red solid, respectively. 
Values and domains of the pattern speeds in solutions at each PA 
are given in Table [3] 

Consequently, solutions with rc=4.1 kpc should yield 
the more accurate description for structure in the zone 
r<4 kpc. This analysis therefore at best indicates that 
within rc=4.1 kpc the data at the nominal PA are well 
fit by two pattern speeds. In addition, though, it seems 
possible to extend the multi-speed model's estimate 
for Slp(r) to 5.3 kpc without loss of validity, and this 
appears to be a good approximation to the pattern 
speeds of the structure across this zone. 

Future high-resolution CARMA observations of M51 
should enable the TWR measurements inside r<4 kpc to 
be more clearly distinguished, especially at the innermost 
radial bins. Presently, however, it is nevertheless clear 
that with the radial calculation at PA=170° we measure 
a pattern speed for the bright spiral structure in the 
zone r<4 kpc higher than the global ^ 38 km s~^kpc~ 
found by iZimmer. Rand fc McGrawl (|2004[ ). also at 
PA=170°. Interestingly, our measurement of a higher 
inner speed resembles the lower bound on such a pattern 
available with the TW calculatio n, f7p>88 km s^^kpc^ 
(jZimmer. Rand fc McGrawll2004f ). 

4.3. Dependence on PA 

We expect the rather large PA uncertainty (5p^=±5° to 
introduce significant variation in the values measured at 
the nominal PA; in the previous section, we used this to 
define the error in the measurement of the best-fit param- 
eters for PA=170°. But TWR solutions from PA=165° 
and 175° themselves additionally indicate a departure 
from the parametrization characteristic of the lowest xt 
solution measured at 170°. 

When we identify the best-fit solutions strictly by their 
x1 over all slices at each iM-rather than restrict our con- 
sideration at PA=165° and 175° to pattern speed solu- 
tions optimal at PA=170°-we find that the values and 
domains of the best-fit pattern speeds vary from PA to 
PA. Figure [S] shows the best-fit solutions at the three 
PAs, the values and xl for which are given in Table [31 
There, solutions with rc=5.3 kpc at PA=165° and 175°, 
unlike at PA=170°, measure at most two distinct pattern 
speeds. More notably, although the best-fit solution with 
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Note. — Representation of the goodness of fit given by x?si ^^ ™ Table [2] for model solutions at 
PA=165° and 175° calculated with rc=4.1 kpc (left) and rc=5.3 kpc (right). Here, the zones 0<|j/|<rt,i, 
'''t,l<\y\<i't,2, and rt,2<|?/|<5.3 kpc are labeled A, B, and C, respectively, with transition radii as identified 
in the lowest-x'^ two (three) pattern speed solution calculated with rc=4.1 kpc (5.3 kpc). At 165°, these 
transitions occur at rt,i=rt=2.3 kpc where rc=4.1 kpc and rt_i=2.3 kpc and rt_2=4.4 kpc where rc=5.3 
kpc. At 175°, rt i=rf=2.3 kpc where rc=4.1 kpc and rt_i=2.3 kpc and rt^2=3.7 kpc where rc=5.3 kpc. 
Pattern speed estimates in units of kms~^kpc~^ in each zone are also listed. For quadratic solutions, the 
values in the first and last radial bins arc indicated. 

best-fit solution with re— 5.3 kpc (4.1 kpc). From the xt s 
it is apparent that, for rc—5.3 kpc two pattern speeds in- 
side r~4 kpc fit the inner two zones significantly better 
than the single pattern speed shown in Figure [H How- 
ever, the xl.s ill the third zone of this triple pattern speed 
solution is rather high; the small 11 kms~^kpc~^ differ- 
ence in the speeds measured inside r~4 kpc therefore 
seems available only at the expense of accuracy in third 
pattern speed. Furthermore, for rc=4.1 kpc the distinc- 
tion between the two pattern speeds measured inside r~4 
kpc weakens; the two in this case are nearly identical to 
the constant value (such that the xt (and xt s) of both 
solutions are comparable). Overall, then, Table [3] sup- 
ports a conclusion that the data at 175° are at best con- 
sistent with a single constant pattern speed inside r~4 
kpc. 

This apparent preference for the measurement of a 
constant pattern speed at PA=175° may suggest a phe- 
nomenological difference in the projection of asymme- 
tries (both intensity and velocity) from that in the 170° 
case. For example, at PA=175°, and also at a higher 
PA=180°, <v> measurements in slices |j/|<2.3 kpc are 
about '^50 % smaller than at PA=170°, a significant dif- 
ference given the flux error for these slices. While this is 
consistent with expectations for the large change in TW 
integrals introduced by a change in the PA (e.g. from 
a nominal SpA=0°; iDebattistal l2003^ ■ the smaller pro- 




FlG. 7.— The best-fit regularized solution for PA=175°-|-/-5° 
and rc=4.1 kpc (5.3 kpc). Dashed red lines represent the difference 
from solutions derived with a one- (two-)pattern speed model at 
PA=170° and 180° (the best fit functional forms at PA=175°). 
Horizontal error bars represent the dispersion in rt calculated in 
solutions with rc=5.3 kpc from PA to PA. Where rc=4.1 kpc, 
the value in the zone of the bright spiral structure corresponds to 
Op^l=59-7/-3 kms~^kpc~ out to rc=4.1 kpc and where rc=5.3 
kpc (shown in thinner line), Hp i=64-2/+15 kms~^kpc~^ out to 
rt=3.7±0.2 and Op,2=ll+2/+9 kms'^kpc"^ out to rc=5.3 kpc. 
Curves for f2, Q±k/2 and Q-k/4 (see ij 14.51 1 are shown in gray. 

rc=4.1 kpc at 165° measures two pattern speeds, at 175° 
no unique pattern speed is measured inside r<2 kpc. 

Model comparisons based on xt over all slices and the 
separated xt.s diagnostic (Table [3]) demonstrate the de- 
gree to which these best-fit solutions differ from those at 
the nominal PA. At 165°, for instance, the xt for solu- 
tions with rc=5.3 kpc suggest that three pattern speeds 
are nearly indistinguishable from the best-fit solution. 
The xt s confirms that the third zone in solutions with 
rc=5.3 kpc is fit equally as well by a third pattern speed 
as by the second speed in Figure [6l (Inside r^A kpc, the 
X^ (and speeds) of the three and two speed solutions are 
nearly identical.) Furthermore, where rc=4.1 kpc two 
pattern speeds fit the data significantly better than a 
single, constant pattern speed. 

At PA=175°, too, judged overall by the xt^ three (two) 
distinct pattern speeds seem to fit nearly as well as the 



jected streaming velocities in this case may be more dif- 
ficult to extract with the TWR calculation, and lacking 
strong signatures, reproducing the higher pattern speed 
measured at 165° and 170° may therefore be improbable 
at 175°. In effect, the single measured speed inside r~4 
kpc at PA=175° may therefore reasonably describe two 
pattern speeds with significant error in each; this speed 
is nearly consistent with the solution plotted in Figure [3] 
with the errors defined by the PA uncertainty. 

Critically, however, the data admit both PA =170° and 
PA=175°, and according to the findings of IShettv et al.l 
(|2007l ) addressed in section 14. 6[ the latter may be ar- 
guably more valid at the inner radii than our chosen 
170°. I n this case, if , as m ight be indicated by the anal- 
ysis of IShettv et al.l ()2007( ). the PA does not reach the 
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assumed 170° until r~3 kpc, rather than measuring a 
distinct pattern speed, fip^i could thus be interpreted as 
simply identifying the region in the disk where the as- 
sumed PA is inappropriate. Furthermore, by the same 
token as above, the large (by comparison) <v> at 170° 
may themselves reflect a misrepresentation of velocities 
in projection from the 175° case. Consequently, if choos- 
ing 170° introduces streaming motions that are unreal, 
the two pattern speeds measured in the best-fit solution 
at 170° could just as persuasively reflect a large PA error 
introduced into the measure of a single constant pattern 
speed. 

Since solutions at the two PAs indicate quite indepen- 
dent radial behaviors, we include here an estimate which 
may be more appropriate for a nominal PA=175°. Fig- 
ure[7]plots solutions with rc=4.1 kpc and rc—5.3 kpc for 
PA=175°-|-/-5° where we have fixed the functional form 
of solutions at PA=170° and 180° to that of the best-fit 
175° solution. (These solutions are also the best-fitting 
at PA=180°; for 170° the multi-pattern speed solutions 
in the previous section are otherwise best). Note that 
the values inside r~4 kpc are only slightly modified with 
the inclusion of the zone 4.1<r<5.3 kpc. 

4.4. The state of current measurements 

Although the PA of the disk is ambiguous, for the par- 
ticular case of a single assumed PA=170°, the majority 
of our analysis leads us to consider the solution in Figure 
[3] a fair representation of the (isolated) inner disk. As 
stated previously, the errors represent PA uncertainty 
introduced to the parameters of the best-fit solution de- 
rived at the nominal PA with a two-pattern speed model. 
This uncertainty 5pa=±5° defines 22% and 16% error on 
the pattern speed estimates flp^i and Op^2, respectively, 
and 14% error in the transition rt, all reasonabl e with 
regard to the standard set by the study of Dcb attistal 
((2003). 

According to the study of iMeidt et~al] ()2008l ). part 
of this error can be expected to have originated with 
limitations in determining the location of the transi- 
tion between the two patterns (assuming they exist), 
as a result of the finite radial bin width. In addition, 
for the inner pattern speed additional uncertainty may 
arise given the disparity between the inner extent of 
the solution and that of the true, dominant two-armed 
pattern, which in the surface brightness terminates at 
the ring-like structure at r'^0.6 kpc. If structure in- 
sid e rcfl kpc, perhaps lik e that identified in the near-IR 
bv lZaritskv. Rix fc Riekel (|1993f ). contributes to the cal- 
culation with a unique pattern speed in this zone, our 
measurement fip.i would represent a combination of this 
value with that for structure out to r~2 kpc. (Note, 
too, in this case, flp^i would also mis-estimate the true 
pattern speed between 0.6<r<2.0 kpc.) Unfortunately, 
identifying whether or not an additional, unique pattern 
exists inside 7'^ 1.0 kpc, or even establishing an inner- 
most extent for the measure flp^i is currently beyond our 
capability; the total degrees of freedom for even the low- 
est order polynomial exceed the number of available bins 
in the innermost radial zone. 

Presently, the pattern speeds in the best-fit solution at 
PA=170° in general tend to be arranged adjacent to the 
angular rotation curve (or perhaps even along the curve 
Q — k/A; see Figure[3]or Figure |4] showing rotation curves 



established in § 14. Sp , much as if identifying a propensity 
towards a material pattern description. Rather than fur- 
nish a description for arms that are material and wind- 
ing, however, we note that the very alignment of multiple 
pattern speeds with the disk angular rotation may re- 
late to an underlying mechanism governing the existence 
and maintenance of structure in the disk. In one inter- 
pretation, the succession of corotation radii implied by 
the best-fit solution might be an indication of resonance 
overlap, as discussed inconclusively in § 14.51 Associated 
with mode-coupling, this would allow quasi-static spiral 
struc ture to be maintai ned over a large portion of the 
disk (ISvgnet et al.lfl988l ) while transporting energy and 
angular momentum outward. 

Our TWR solutions furthermore seem unlike what 
might be expected for transient density waves, with de- 
scription deriving f rom the propagation of tid al pertur- 
bations studied by ISalo fc LaurikainenI l2000bl For ex- 
ample, although their range of applicability seems lim- 
ited to the innermost radii (but taken as an approxima- 
tion to the best-fit multi-speed solutions) the bin val- 
ues in the current set of quadratic solutions are much 
closer to the angular rotation of the disk than 17 — k/2, 
near which much of the m= 2 structure in the models of 
ISalo fc LaurikainenI (|2000bD achieves its greatest ampli- 
tude. 

In order to best establish the extent to which the TWR 
solution in Figure [3] is truly a valid description of the 
bright spiral structure, in the section immediately follow- 
ing, and in §§ l4.6l to l4.71 we relate the radial dependence 
exhibited by the solution to observed morphological and 
kinematic structure. The inner disk of M51 has been sug- 
geste d to sustain radial variation in t he PA (|Shetty et"an 
I2007t ) and an additional m=3 mode (jHenrv et al.ll2003f l. 
both of which undeniably challenge the authenticity of 
the TWR solutions, and so we address the possible in- 
fluence of each of these in turn. 

4.5. Possible complimentary evidence for multiple 
pattern speeds and indications of mode coupling 

Though perhaps unexpected, the identification of a 
transition between two pattern speeds in the inner disk 
seems supported by independent studies of the bright spi- 
ral structure. At least two sections best fit with slightly 
different pitch angles ip have been identified in both spi- 
ral arms, possibly the signature of two or more distinct 
pattern s peeds. Notably , the anisotropic wavelet ap- 
proach of iPatrikeev et all ()2006l ) shows evidence for ex- 
treme departures from the conventionally adopted value 
ip=21°. The maximum occurring nearly symmetrically 
in both arms at r~2 kpc (see Figures 6, 7a and 8a in 
IPatrikeev et al.l f2006f )-verv near the transition rt iden- 
tified in our best-fit solution at PA=170°-is especially 
compelling since it may indicate more than a simple de- 
parture from a logarithmic dependence. The transition 
rt.2 in Figure[3]also occurs near a maximum in ip. (To be 
sure, the other extrema in ip imply no such correlation) . 
This analysi s is la rgely consistent with the findings of 
I Henry et al.l (|2003[ ) covering radii r<4 kpc which iden- 
tify three arm sections each with a unique ip. 

Of course, a systematic error in deprojection, such 
as due to an incorrectly adopted PA or inclination an- 
gle (or a radially varying PA, as considered in § 14. 6p 
could very well alone produce the effect measured by 
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iPatrikeev et all (|2006f) (where PA=170°). A firm con- 
clusion might therefore require a better understanding 
of how a change in pitch angle at a given radius relates 
to a change in pattern speed, for instance (assuming that 
the spirals are indeed logarithmic). 

The transition between two patterns inside r ~4 kpc 
indicated in our solution also seems significant given that 
it coincides with features in the zeroth moment map's 
Fourier decomposition (Figure [2]); as at r^A kpc, the 
power in the m=2 mode is characterized by a decline at 
r~2 kpc possibly marking the termination of a distinct 
structure. (The same can be inferred at 4 kpc which co- 
incides with the transition rt_2 in the three-pattern speed 
solution.) 

Perhaps more compellingly, the transitions parameter- 
ized in our solutions at 170° also appear to coincide with 
resonances, as illustrated in Figures [3] and [H This, of 
course, would seem to depend largely on the assume d 
rotation curve. As first demonstrated bv iTuUvl (|1974l ). 
streaming motions appear non-negligibly in the rotation 
curve of M51, making the true circular vel ocity difficult 
to con strain. As more recently cataloged bv lShettv et al.l 
(|2007t l. all other kinematic parameters are likewise sus- 
ceptible to such errors, and so rotation curves gener- 
ated with them are prone to inaccuracy. To reduce the 
impact of streaming motions (and perhaps other sys- 
tematic errors) on our resonance identifications, we fit 
our own ROTCUR-derived rotation curve with the com- 
monl y used approximation (e.g. by iFaber fc Gallagherl 
[l979h 



Vrot{r) 



''max y I '^max } 



{l/3 + 2/3{r/r,na.rf^" 



(5) 



which yields a smoothed curve for Q. In this expression, 
Vmax is the maximum rotational velocity, r^ax is the lo- 
cation where Vmax occurs, and n determines how rapidly 
the curve becomes Keple rian. Alternative fits (like that 
used on the inner 30" bv lAalto et al.|[l999l ) supply simi- 
lar conclusions. 

With the resulting curves for 17, fi — k/2, fl — k/4, 
and n + k/2 plotted in Figures [3] and [4] we highlight 
the possible locations for the corotation, inner Lindblad, 
inner 4:1 ultraharmonic, and outer Lindblad resonances 
(or CR, ILR, UHR, and OLR) for each measured pat- 
tern speed. Immediately we notice that both pattern 
speeds in the solution with rc=4.1 kpc end at their CR 
within the uncertainties. This circumstance is consis- 
tent with an e arly prediction for where spirals terminate 
(i.e. lLinlll97Cil ). which later yielded to findings that spi- 
rals can e xtend as far as OLR, if sometimes faintly (see 
[Elmegrecn et al.lll989l for example). 

In addition, the transition between the two pattern 
speeds appears to occur at a resonance overlap. As 
demonstrated in Figure [3l the CR of fipi overlaps 
the UHR of the pattern with flp^2- Such coincidences 
have been identifi ed in barr ed spiral s imulatio ns o f 
iRautiainen fc Said (|1999D and iDebattista et all (|2006l ). 
As the former investigate, this overlap at resonance 
may be characteris tic of non-linear mode coupling (e.g. 
iTaggeret aLlll987l and lSvgnet et al.lll988[ ) whereby en- 
ergy and angular momentum are transferred between the 
mo des. But in contr a st to the CR-ILR overlaps studied 
by iMasset fc Tagged ()1997t ) , which are accompanied by 
boosted beat modes detectable in the simulation power 



spectra at the overlap, they find no comparable evidence 
for mode-coupling in the case of the CR-UHR overlap. 
(They suggest this overlap may nevertheless be related 
to a physical process.) 

A CR-ILR overlap between ilp^i and Qp^s in the so- 
lution with rc=5.3 kpc, on the otherhand, may be vi- 
able within the uncertainties. However, between flp^2 
and f2p^3 in Figure |4] a similar resonance overlap is not 
so clear; near the transition r^ 2 CR of f2p 2 falls between 
the ILR and the UHR of f^p^g! 

Figure [3] also exhibits a turnover in the curve f2 — k/2, 
suggesting that patterns with angular speeds above the 
maximum lack an ILR. However, given uncertainty in the 
rotation curve, this is diffic ult to constr ain: while the 
angular frequency curves of iTullvl (|1974f ) indicate that 
fl — k/2'^4 7 kms~^kpc~ at the turno ver (also repro- 
duced in the lSalo fc LaurikainenI (|2000bl ) model), we find 
that the maximum occurs at ~75 km s~"'^kpc~ . (Our fit 
for ri(r) may be slightly steep inside ~1.0 kpc.) Never- 
theless, it is apparent for this solution that flp^i lacks an 
ILR. This suggests that a (trailing) wave with fJp^i can 
reflect from the center as a l eading structure, a ci rcum- 
stance complimentary to the iScoville et al.l (|2001l ) HST 
observations of central lead ing waves, as pointed out by 
ISalo fc LaurikainenI (|2000bD . 

4.6. Effect of a radial variation in PA 

From their analysis of the CO and H^ kinematics, 
iShettv et al.l ([^007) find evidence for large non-zero ra- 
dial flux (as measured by the mass/surface brightness- 
weighted radial velocity) in radial ranges that depend on 
the c hoice of PA. From F igure 18 of that study, in partic- 
ular, IShettv et al.l (|2007D speculate that mass flux could 
be conserved should the disk of M51 sustain a radially- 
dependent PA (and/or inclination). This could be ap- 
proximately achieved with PA=180° out to r'^1.8 kpc, 
PA=175° out to r-2.8 kpc and PA=170° out to r-3.7 
kpc. (The inclination angle, which might also be ex- 
pected to vary, is much harder to account for in the TWR 
calculation.) 

If the PA does vary radially then the measurement in 
Figure |3| (or Figure [7]) could be affected by projection 
errors at certain radii. Note that a radially varying PA 
implies a warp in the inner disk (which, if real, could 
be due to the presence of the companion) and so the 
disk would also not meet the assumptions of the TWR 
method. Interpreted in this manner, our finding of pos- 
sible multiple pattern speeds in the inner disk may be 
the result of such an effect. 

We explore this possibility by allowing the PA to vary 
radially in the TWR quadrature according to the pre- 
scription given at the beginning of this subsection. For 
simplicity, we retain i—2A° throughout the disk and let 
PA=170° at all radii beyond r'^3.7 kpc. Figure |5| shows 
the best-fit solutions with rc=4.1 kpc. Errors represent 
a residua l PA uncertainty 6p A—^2° estimated from Fig- 
ure 18 of lShettv et"an (|2007f ). 

Interestingly, the global pattern speed inside rc=4.1 
kpc (r2^^=62±2 kms~^kpc~ ) closely resembles the 
measurement at PA=175° (Figure [7]). The best-fit so- 
lution with rc=5.3 kpc (0^;[^=62±2 km s~-^kpc~"^ out to 
rt=3.8±0.5 kpc, f7^^=18±2 km s'^kpc"^ out to rj also 
resembles that at 175°, and here the pattern speed in the 
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Fig. 8. — The best- fit regularized solution for the PA twist map 
with 5pj^=±2° and rc=4.1 kpc (5.3 kpc) shown in thick (thin) 
line. (See text for a description.) Dashed red lines represent the 
difference from solutions derived with a two-pattern speed model 
best for the SpA=0° case; errors for the solution with rc=4.1 kpc 
are nearly coincident with these lines and have been left off for 
clarity. Horizontal error bars represent the dispersion in rt calcu- 
lated in solutions with rc=5.3 kpc. Where rc=4.1 kpc, the value in 
the zone of the bright spiral structure corresponds to f!^ =61-1-4/- 

1km s~^kpc~^ out to rc=4.1 kpc and where rc=5.3 kpc (shown in 
thinner line), f2^^=62it2 km s~^kpc~^ out to rt=3.8it0.5 kpc and 
Slp2=18i2 km s~ kpc~ out to rc=5.3 kpc. For there solutions, 
bins exterior to re (not shown) are calculated without regulariza- 
tion. 

zone 4<r<5 kpc does not seem to be the result of an in- 
correctly assumed PA (i.e. PA=170° instead of 165°); 
extending the twist by another 5° at radii greater than 
r—3.7 kpc produces httle change in the calculated solu- 
tions. 

That the estimates in Figures [S] and [7] are so similar 
seems to suggest the twist solution is less a manifestation 
of the PA twist than an indication that the PA assumed 
here inside r~3 kpc is everywhere closer to 175° than 
170°. The mean PA of the twist is 175°, so it may be 
reasonable to infer that the twist solution predominantly 
reflects information nearly identically to the 175° case. 
Solutions at 180°, too, are very similar to those at 175°, 
as indicated by the estimate assembled in Figure [71 

Since imposing the twist does not seem to introduce a 
novel character to the TWR measurement, by extension 
this leads us to conclude that the regularized TWR cal- 
culation is insensitive to minor radial variation in the PA 
(i.e. (5pyi=±5° over roughly 4 kpc), if real. However, if 
the disk PA is assertably closer to 175° than 170°, this 
seems to reinforce the impression that the bright spiral 
structure may be best described by a single constant pat- 
tern speed. 

Though compelling, we emphasize that this exercise 
should not be interpreted as confirmation or denial of ra- 
dial variation in the PA, nor as providing an unequivocal 
measure for the pattern speed of the bright spiral struc- 
ture. Critically, imposing the twist tends to remove the 
most noticeable asymmetries in the velocity field, par- 
ticularly within 60" , very much in the manner described 
previously for PA=175°. If true signatures of pattern 
speeds have been obscured or eliminated at this PA, this 
may prevent the measurement of a distinct pattern speed 
inside r^2.0 kpc. 

4.7. Relation to m=3 structure 



In principle, TWR measurements at either PA=170° or 
175° may reflect signatures of patterns other than those 
of the bright two-armed structure alone. The Fourier 
power spectrum of the surface brightness reveals rich 
structure in the disk of M51, much of it coexisting over 
roughly 2 kpc in the inner disk. In this section we review 
the part icular possibi l ity th at weak m=3 structure iden- 
tified bv lRix fc Riekd (|1993D out to ^^2.5 kpc contributes 
with a unique pattern speed to out TWR solutions. 

Only if all structures in the same radial zone have iden- 
tical pattern speeds will the TWR calculation accurately 
reflect this sole speed; given a surface brightness distri- 
bution which reflects two coincident contributions (say, 
from 771=2 and m=3 structures) with unique time depen- 
dence (i.e. different pattern speeds), the TWR calcula- 
tion is currently unequipped to constrain either one or 
the other. Though a generalization can be made under 
the assumption that both pattern speeds are constant, 
it is beyond the scope of this work to develop either the 
TW or TWR calculation appropriate to the situation. 

For M51, it may be possible in the future to directly 
relate the velocity asymmetry from arm t o arm of the 
bright m=2 spiral structure identified by iHenrv et al.l 
()2003f ) to the presence of the m=3 mode, an d to its pat- 
tern speed in particular. (|Henrv et al.ll2003l have already 
successfully demonstrated that the presence of the to=3 
mode induces a systematic offset in the azimuthal po- 
sitions of the two main arms). This should allow us to 
establish the expected combination of speeds in the TWR 
calculation; if the implied m—3 pattern speed is unlike 
the measure ftp^i found with PA=:170°, for example, this 
speed is presumably unshared by the r7i=2 structure. 

5. CONCLUSION 

In this paper we present regularized TWR solutions 
for the pattern speed of the bright spiral structure in the 
inner disk of M51, derived with velocity and intensity 
information from the ISM-dominant molecular compo- 
nent traced by high-resolution CO observations. These 
solutions are arrived at by isolating the inner disk from 
errors which evidently originate with both the quality of 
sampling/detection and the pattern speed-modeling in 
bins covering the outer, material pattern. So although 
our procedure prevents us from constraining the outer- 
most pattern speed, calculating the outer bins without 
regularization in principle improves the accuracy with 
which the solution for the inner disk can realize the true 
pattern speed. 

Our primary result with this implementation is 
the measurement inside 4 kpc of two pattern speeds, 
both significantly higher, and together fitting the 
data better, than the c o nstant global measure of 
IZimmer. Rand fc McGrawl ()2004[ ) at the nominal 
PA=170°. A third, lower pattern speed, extending 
beyond 4 kpc out to (at least) 5.3 kpc is also detected, 
nearer the speed expected for the material pattern. 
Significantly, the transitions between the measured 
pattern speeds coincide with resonances; the two pattern 
speeds inside 4 kpc both end at corotation within the 
uncertainties. Since it is in no way imposed by the 
method, this dynamically reasonable scenario tends to 
give us confidence as to the physical plausibility of the 
pattern speeds returned by the analysis. 

Of course, given that a pattern speed interior to r~2 
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kpc is only weakly detected (if at all) at PA=175°, the 
accuracy of the description provided by two pattern 
speeds may depend on whether PA=170° or PA=175° 
is mor e accurate, an uncertainty raised by IShettv et al.1 
()2007t ). for instance. If the disk is best described with 
PA=175°, we find evidence that a single constant 
pattern speed inside 4 kpc best characterizes the bright 
spiral structure. Furthermore, as contemplated in § 
4.61 a radially varying PA which decreases from 180° 
(and reaches 170° near 3 kpc), perhaps suggested by 
the results of IShettv et al.l (|2007t ). also favors a single 
measured pattern speed interior to 4 kpc. 

Again, however, while the analysis presented here 
cannot resolve the question as to which PA is more 
appropriate, we find meaningful, independent evidence 
in favor of the pattern speeds measured at 170°, in 
particular. For example, consistent with expectations 
of leading struc ture at the inner m ost radii (as in the 
observations of IScoville et al.l 120011 ). the higher speed 
inside r~2 kpc lacks an ILR. In addition, attendant to 
our finding that both speeds interior to 4 kpc terminate 
at corotation, the transition between the two roughly 
coincides with an inferred location of resonance overlap 
wherein the inner's corotation resonance and the outer's 
inner 4:1 resonance align. The radial domain of the 
pattern speed measured at PA=175°, in contrast, is 
not as clearly associated with resonance radii. Since 
the bright spiral structure does not appear along the 
minor axis near ~2 kpc, the corotation resonance at this 
location implied by the solution at 170° is unfortunately 
unconfirmable through inspection of radial streaming 
velocities under the density wave interpretation. 

We also find remarkable agreement between the char- 
acteristics of the two speeds inside 4 kpc at PA=170° 
and other evidence in the inner disk consistent with 
multiple pattern speeds. The transition parameterized 
in our best-fit solution for the inner disk coincides with 
significant variation in the two-armed spiral pitch angle. 
Since a change in the pitch angle is expected to be 
accompanied by a change in streaming motions, both 
parameters are presumably attendant to the signatures 
(streaming or otherwise) of the patterns. 

Although the pattern speed interior to 2 kpc in the 
solution at 170° (or 175°) may reflect a u nique contri- 
bution from the m=3 mode observed by iRix k, Riekd 
(|1993f ). as described in § 14. 7[ the measurements in 
Figure |4] (or Figure [7]) presumably directly relate to the 
patterns present in the disk and so (depending on the 
PA) should provide a fair description of the dynamics 



therein. As such, it may be possible that ensuing 
observations and studies better discriminate between 
the two seemingly disparate radial dependencies implied 
for the PAs considered here. 

Even at present our TWR solutions yield interpreta- 
tions with which to observationally address the question 
of spiral longevity. That the solutions at both 170° 
and 175° feature constant pattern speeds would imply 
that our solutions are indicative of long-lasting spiral 
structure. Interestingly, at the in nermost radii both 
qualita tively resemble the model of ISalo fc LaurikainenI 
(|2000bD where, characteristic of the isolated evolution of 
the disk, the dominant m=2 component has a constant 
pattern speed ~50 kms~^kpc~ out to ^-^1.2-1.8 kpc. 
As for the region between ~1.8 kpc and 4.6 kpc in those 
models where interaction with the companion introduces 
a succession of transient structures, our solutions at 
both 170° and 175° otherwise describe at most two 
steady patterns in distinct radial zones. 

In the immediate future, observations with higher 
resolution and sensitivity should afford TWR calcu- 
lations with finer radial bins, thereby allowing for 
the parameterization of more distinct radial zones, if 
present. This will either confirm our solutions for rip(r) 
or perhaps demonstrate that solutions describe a succes- 
sion of many disc rete patterns (simila r to the transient 
structures in the ISalo fc LaurikainenI (|2000bO models), 
or simply a winding, material pattern. Again, though, 
our multiple-speed and other, quadratic solutions in 
general more closely follow J7 throughout the disk than 
$7 — k/2 characteristic of m= 2 structure in the models 
of ISalo fc LaurikainenI (|2000bf) . 

Despite the lingering ambiguity in the PA, these 
TWR solutions present a new picture of the bright 
spiral structure of M51, one that should prompt tests 
of long-lived density wave theories in other nearby 
grand-design spirals. At the very least, this study marks 
a successful starting point for continued tests of the 
relation between multiple spiral pattern speeds in a 
single disk; investigations into the number and radial 
domains of pattern speeds and spiral winding in nearby 
spirals will be the subject of upcoming work. 
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